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We extend the classical Landau theory for rotating nuclei and show that the backbending in 162 Yb, 
that comes about as a result of the two-quasiparticle alignment, is identified with the second order 
phase transition. We found that the backbending in 156 Dy, caused by the instability of 7-vibrations 
in the rotating frame, corresponds to the first order phase transition. 
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Quantum phase transitions (QPTs), that occur at 
zero temperature as a function of some nonthermal con- 
trol parameter, are studied extensively in diverse areas 
of physics including, for example, low-dimensional sys- 
tems in condensed matter physics atomic nuclei and 
molecules Q. In nuclear structure, QPTs associated with 
transitions between different shapes or geometric con- 
figurations of a chain of nuclei in the ground state are 
being studied intensively during the past years within 
the interacting boson model (cf Q). Thanks to novel 
experimental detectors, a new frontier of discrete-line 7- 
spectroscopy at very high spins has been opened in the 
rare-earth nuclei. These nuclei can accommodate the 
highest values of the angular momentum, providing one 
with various nuclear structure phenomena. Evidently, 
from analysis of the rotational states one expects to ob- 
tain further insight into the nature of QPTs in finite 
quantum systems. 

One of the well known phenomena at high spins is a 
backbending. A sudden increase of a nuclear kinemati- 
cal moment of inertia ^ 1%> — I/O of the lowest (yrast) 
level sequence as a function of a rotational frequency O 
(see experimental results for 156 Dy and 1G2 Yb in Figfl]) 
is a paradigm of structural changes in a nucleus under 
rotation. While one observes a similar picture for the 
backbending in the considered nuclei (see FigUJi,b), a 
different response of a nuclear field upon the rotation 
becomes more evident with the aid of the experimental 
dynamical moment of inertia Q^ 2 ' = dl/d£l A/AE^ 
as a function of the angular frequency (see Fig[TJ c, d). 
Indeed, the dynamical moment of inertia, due to the re- 
lation Sj( 2 ' = SJW + Qd'<s^ /dfl, is very sensitive to struc- 
tural changes of a nuclear field. At the transition point 
3^ 2 ) wildly fluctuates with a huge amplitude in 156 Dy, 
whereas these fluctuation are quite mild in 162 Yb. One 
of the purposes of the present contribution is to inter- 
pret the experimental results in 156 Dy and 162 Yb (see 
FigfT]) from perspective of the Landau theory of phase 
transitions within a microscopic approach [J] based on 
the cranked Nilsson model plus random phase approx- 
imation (CRPA). We will show that the interplay be- 



tween the alignment (single-particle degrees of freedom) 
and collective quantum fluctuations determines the type 
of the shape-phase transition at the backbending. 

The phase transition is usually detected by means of 
an order parameter as a function of a control parameter. 
In rotating nuclei one can suggest a few order parameters 
like deformation parameters of a nuclear effective poten- 
tial, P and 7, that characterize the geometrical configu- 
ration (cf 0), - as a function of the rotational frequency, 
i.e., the control parameter. To analyze the experimental 
data of above nuclei, we use a cranked Hamiltonian 



Ha = H — 
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The term H = H^n + -ff ac id contains the Nilsson Hamil- 
tonian Hmu and the additional term that restores the 
local Galilean invariance of the Nilsson potential in the 
rotating frame. The Nilsson potential naturally incor- 
porates the quadrupole deformation parameters (/? and 
7) of a nuclear shape. The interaction includes sep- 
arable monopole pairing, double stretched quadrupolc- 
quadrupole (QQ) and monopole-monopole terms. The 
details about the model Hamiltonian ((TJ can be found in 
[U, [E] ■ I n °ur approach mean field parameters are deter- 
mined from the energy-minimization procedure (see 0]). 
The consistency between the mean field and the resid- 
ual interactions of the Hamiltonian (UJ) was achieved by 
varying the strength constants of the pairing and QQ 
interactions in the CRPA. It results in the separation 
collective excitations from those that are related to the 
symmetries broken by the mean field. Among them are 
the conservation of particle numbers and space symme- 
tries (see details in 

We found that the triaxiality of the mean field sets in 
at the rotational frequency Ml c = 0.250, 0.301 MeV for 
162 Yb, 156 Dy, respectively, which triggers a backbending 
in the considered nuclei (see Fig.l in [4J). To elucidate 
the different character of the shape transition from axi- 
ally symmetric to the triaxial shape and its relation to 
a phase transition, we consider potential landscape sec- 
tions in the vicinity of the shape transition. Since we 
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FIG. 1: Rotational behavior of the experimental, kinematical = I/Q and dynamical Si' 2 - 1 ~ 4/A_B 7 moments of inertia. 
Here, hQ = E-,/2, i? 7 is the 7-transition energy between two neighboring states that differ on two units of the angular 
momentum I and Ai? 7 is the difference between two consecutive 7-transitions. The experimental data denoted by black 
squares are taken from ||. The experimental rotational frequency at the transition point is hQ c ~ 0.27,0.32 MeV for 162 Yb 
and 156 Dy, respectively. The results of calculations for Sr 1 ' are connected by a solid line. 



analyze a shape transition from the axially symmetric 
shape (7 = 0) to the triaxial one (7 7^ 0), we choose only 
the deformation parameter 7 as the order parameter that 
reflects the broken axial symmetry. Such a choice is well 
justified, since the deformation parameter (3 preserves its 
value before and after the shape transition in both nuclei: 
Pt « 0.2 for 162 Yb and (3 t « 0.31 for 156 Dy. Thus, we 
consider a mean field value of the cranking Hamiltonian, 
En(l',Pt) = (Hn), for different values of SI (our state 
variable) and 7 (order parameter) at fixed value of (3f 

For 156 Dy we observe the emergence of the order pa- 
rameter 7 above the critical value Ml c = 0.301 MeV of 
the control parameter SI (see a top panel in Fig[2]). Below 
and above the transition point there is a unique phase 
whose properties are continuously connected to one of 
the coexistent phases at the transition point. The order 
parameter changes discontinuously as the nucleus passes 
through the critical point from axially symmetric shape 
to the triaxial one. The polynomial fit of the potential 
landscape section at Ml c = 0.301 MeV yields the follow- 
ing expression 

F(0; 7 ) = F (ft) + F 2 (ft)7 2 -F 3 (fth 3 + ^4(r>)7 4 , (2) 

where the coefficients Pb(0) = 0.3169 MeV, 7 in de- 
grees and F 2 (n) = 0.12239, F 3 (fi) = 0.009199, F 4 (fi) = 
1.7 x 10 -4 are defined in corresponding units. We can 
transform this polynomial to the form 



F(n;<y)-F (Sl) rf 
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(3) 



where F = (3F 3 ) 4 /(4F 4 ) 3 , a = 8F 2 F 4 /(9F 3 2 ), r, = 
4F 4 /(3F 3 )7. The expression ^ represents the generic 



form of the anharmonic model of the structural first order 
phase transitions in condensed matter physics. The con- 
dition dF/drj = determines the following solutions for 
the order parameter r\ : 7^ = 0, r\ 2 ,z = (1± \/l — 4a)/2. 
If a > 1/4, the functional F has a single minimum at 
7/ = 0. Depending on values of a, defined in the interval 
< a < 1/4, the functional F manifests the transition 
from one stable minimum at zero order parameter via 
one minimum+metastable state to the other stable min- 
imum with the nonzero order parameter. In particular, 
at the universal value of a = 2/9 the functional F has 
two minimum values with _F = 0atr7 = 0— > 7 « 0° 
and 77 = 2/3^7« 27° and a maximum at r\ — 1/3 — > 
7 13.5°. The correspondence between the actual value 
7 « 20° and the one obtained from the generic model is 
quite good. 

In general, the alignment of angular momenta of a nu- 
cleon pair occupying a high-j intruder orbital near the 
Fermi surface is considered as a main driving force that 
leads to the backbending. Although two-quasiparticle 
states align their angular momenta along the axis x (col- 
lective rotation), the axial symmetry persists till the 
transition point (see also discussion in Q). The RPA 
analysis of 7-vibrational (lowest) excitations of the pos- 
itive signature in the vicinity of the shape transition 
demonstrates a collective nature of these excitations. The 
mode blocks a transition to the triaxial shape. However, 
at the transition point, this mode is anomalously low in 
the rotating frame. It appears that soft positive signa- 
ture 7-vibrations (fluctuations of the order parameter 7) 
coupled to the other modes are responsible for the shape- 
phase transition of the first order. A drastic change of 
the mean field configuration leads to large fluctuations of 
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FIG. 2: The rotational dependence of the order parameter 7 and the energy surfaces sections F(Q, 7) = Eni'y; fit) — E m in for 
156 Dy (top) and 162 Y b (bottom) before and after the transition point. The energy is given relative to the value E m i„ = En{/3t,y) 



at Ml = 0.255, 0.302 MeV for 162 Yb, lba Dy, respectively. 
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the dynamical moment of inertia at the transition point, 
since 3^ 2 ^ = —d 2 Efi/d 2 £l, which is reproduced success- 
fully in the CRPA (see details in Q). It seems reasonable 
to say that the backbending in 156 Dy possesses typical 
features of the first order phase transition. 

In the case of 162 Yb the energy Ea(j; f3 t ) and the 
order parameter (Fig(2J) are smooth functions in the 
vicinity of the transition point f2 c . This implies that 
two phases, 7 = and 7 7^ 0, on either side of the 
transition point should coincide. Therefore, for fl near 
the transition point O c we can expand our functional 
F(Q,j) = J5n(7)A) - E mm in the form 

F(fi; 7) = Fi(fi)7 + ^(^h 2 + F 3 (Q)j 3 + F 4 {n)j 4 + ... 

(4) 

The energy surfaces are symmetric with regard of the 
sign of 7 and this also supports the idea that the effec- 
tive energy F can be expressed as an analytic function of 
the order parameter 7. With the aid of the conditions of 
the phase equilibrium, dF/d-f — and d 2 F/d-f 2 > 0, one 
can show that Fi(il c ) — i*2(^c) = 0. In virtue of these 
equations and of the fact that all phases at the transi- 
tion point should coincide, we obtain from dF/d"f = 
that F 3 (Q = Q c ) = 0. Assuming that F 3 = for all SI, 
the minimum condition dF/dj — yields the following 
solution for the order parameter 

F 2 (0) _ f ^0 foril^n 



71 = , 72,3 



2F 4 (fi) 



= for il = Vl c 



(5) 



Since at the transition point F2(£l c ) = 0, one can propose 
the following definition of the function ^(fi): 



(6) 



Thus, we have 7 ~ (f2 — £l c ) L ' and the critical exponent 
v = 1/2, in accord with the classical Landau theory, 



where the temperature is replaced by the rotational fre- 
quency. Indeed, the numerical results (see Fig|2]) are in 
an agreement with Eqs.([5]),([n]): 7 = for Ml < Ml c , 
while 7 7^ for Ml > Ml c . Thus, the backbending in 
162 Yb can be classified as the phase transition of the sec- 
ond order. The smooth behavior of the function F at 
the transition point implies a small amplitude of fluc- 
tuations of the dynamical moment of inertia. This re- 
sult is nicely reproduced within the CRPA approach 
Q. The RPA analysis of the lowest 7- vibrational mode 
in 162 Yb indicates on the breakdown of the quadrupole 
phonon. At the vicinity of the transition point one proton 
and one neutron two-quasiparticle components dominate 
(~ 95%) in the phonon structure and the backbending is 
caused by the alignment of the neutron two-quasiparticle 
configuration. 

The identification of wobbling (negative signature) ex- 
citations near the yrast line provides sure evidence of the 
onset of the triaxiality. We found (see [H) that a shape- 
phase transition produces relatively high-lying wobbling 
vibrational states in 156 Dy. In contrast, a soft shape- 
phase transition from the axially deformed to nonaxial 
shapes produces the low-lying wobbling excitations in 
162 Yb. We have also established the relation between 
the sign of the 7-deformation and selection rules for the 
quadrupole electric transitions from the wobbling to the 
yrast states (see Q). From our calculations it follows 
that at low angular momenta {Ml < 0.28, 0.3 MeV in 
162 Yb, 156 Dy, respectively) the first negative signature 
one-phonon band populates with approximately equal 
probabilities the yrast states with I' = I ± 1 (I is the 
angular momentum of the excited state) . At Ml c a shape- 
phase transition occurs, that leads to the triaxial shapes 
with the negative 7-deformation in the both nuclei. In 
turn, the negative signature phonon band decays stronger 
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on the yrast states with angular momenta I' = I — 1, 
starting from Ml > M7 C . We predict also the dominance 
of AJ = 1H magnetic transitions from the wobbling to 
the yrast states, independently from the sign of the 7- 
deformation (see |S|). 
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